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Abstract. We study the mechanics of a reversible decohesion (unzipping) of an elastic layer 
subjected to quasi-static end-point loading. At the micro level the system is simulated by an 
elastic chain of particles interacting with a rigid foundation through breakable springs. Such sys- 
tem can be viewed as prototypical for the description of a wide range of phenomena from peeling 
of polymeric tapes, to rolling of cells, working of gecko's fibrillar structures and denaturation 
of DNA. We construct a rigorous continuum limit of the discrete model which captures both 
stable and metastable configurations and present a detailed parametric study of the interplay be- 
tween elastic and cohesive interactions. We show that the model reproduces the experimentally 
observed abrupt transition from an incremental evolution of the adhesion front to a sudden com- 
plete decohesion of a macroscopic segment of the adhesion layer. As the microscopic parameters 
vary the macroscopic response changes from quasi-ductile to quasi-brittle, with corresponding 
decrease in the size of the adhesion hysteresis. At the micro-scale this corresponds to a transition 
from a 'localized' to a 'diffuse' structure of the decohesion front (domain wall). We obtain an 
explicit expression for the critical debonding threshold in the limit when the internal length 
scales are much smaller than the size of the system. The achieved parametric control of the 
microscopic mechanism can be used in the design of new biological inspired adhesion devices 
and machines. 

Key words: Unzipping, Adhesion, Peeling, Hysteresis, DNA, Gecko, Calculus of Variations, V- 
convergence. 



Introduction 

Adhesion phenomena are governed by complex energy exchanges between multiple scales 
representing hierarchical structures. The phenomenological modeling of adhesion has been suc- 
cessful in describing a variety of experimentally observed static and dynamic regimes (see 
[9j [131 123"! d]). The phenomenological models, however, are of a black box type and have lim- 
ited predictive power outside of a particular range of parameters. More importantly, they can 
not be used for the microstructural optimization of the artificially created adhering materials 
and nanorobotics devices. It is therefore of interest to develop a multi-scale approach linking 
the microscopic attachment-detachment mechanisms with the macroscopic phenomenological 
parameters. This step is crucial for the analysis of a variety of adhesion related phenomena from 
fiber decohesion in composites [T7] and crazing/peeling phenomena in polymers [27\ 113] . 
to the activity of focal adhesions involved in cell motility [8j [181 E] an d the low temperature 
denaturation of single molecule DNA [22J. 
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In this paper we contribute to this general task by considering the minimal model ac- 
counting for the interplay between elasticity and cohesion. While this model captures only the 
most important effects associated with quasi-static decohesion, it has the advantage of being 
amenable to a completely transparent mathematical analysis allowing one to study how macro- 
scopic responses vary depending on the microscopic parameters. We focus on the case when the 
debonding is reversible. This type of decohesion, also known as unzipping, is crucial for the work- 
ing of a variety of biological systems (|_26j), in particular, for the functioning of the self-similar 
fibrillar Gecko's hair ([12 ^ 121 ^ 1301). The irreversible version of the model, which is more adequate 
for the description of decohesion in conventional engineering materials, will be presented in a 
separate paper. For a general comparison of reversible and irreversible fracture see [6]. 




FIGURE 1 . Finite elastic layer attached to a rigid background. The controlling parameter 
is the displacement d applied at one end, the other end is free. The measured response is 
represented by the dependence of the point force A on d. 

According to typical observations the process of reversible decohesion due to quasi-static 
point loading includes three main stages ( j!3j and references therein). First the system behaves 
elastically until decohesion begins. Then there is a steady state incremental evolution of the 
decohesion front. Finally, at a critical threshold, the system undergoes a sudden transition to 
the fully debonded configuration. Similarly, if the system is unloaded from the fully debonded 
state, there exists an unloading threshold beyond which a finite portion of the adhesive layer 
suddenly reattaches to the adhesive background. The whole phenomenon is usually hysteretic 
with different attachment and detachment thresholds. 

In an attempt to reproduce this basic behavior we consider a chain of massless points con- 
nected by harmonic shear springs. The particles are attached to a rigid support by breakable 
elastic links (see FigQ]) which mimic, depending on the parameters, either direct molecular inter- 
actions or interactions through the fibrillar adhesive layer with internal elasticity. For simplicity 
we neglect the bending stiffness of the elastic layer which could be accounted for by adding elas- 
tic interaction of next to nearest neighbors [16} [20] . We consider a quasi-static transversal point 
loading of the otherwise free layer in a hard device and study the rate independent evolution 
of the emerging debonding front which can be viewed as a domain wall separating bonded and 
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debonded phases [22J. The reversal of the front direction represent the switch between zipping 
and unzipping. 

Similar discrete models of the Frenkel-Kontorova type have been used previously in the 
analysis of lattice trapping of cracks in crystals [14], interfacial wetting [15] and other phenomena 
where an on site potential with sublinear growth competes with an elastic coupling of the 
elements [2] . In connection to the duplication and transcription of the DNA an approach of this 
type was first proposed by Peyrard and Bishop who applied it to the modeling of equilibrium 
melting phase transition (see the review [22]). Our model can be viewed as a purely mechanical 
version of the Peyrard-Bishop model where we go beyond global minima of the energy (zero 
temperature limit) and investigate the structure of the whole energy landscape (see also |28j). 
Our use of the simplified piece-wise quadratic approximation for the on site potential allows us 
to formulate the main results in the analytic form. 

We use an incremental energy minimization approach and solve the finite dimensional 
variational problem for each value of applied displacement. Due to the simplicity of the cohesive 
potential we are able to find all equilibrium configurations and identify those representing global 
and local minima of the energy. Our analysis shows that the local minimizers of the energy 
always have a single decohesion front. The metastable configurations, forming separate branches 
parameterized by the loading parameter, can abruptly end. The absence of continuity leads to 
the necessity of 'dynamic snapping' from one branch to another. While these events may be 
dissipative, they do not prevent overall reversibility (see also \12 \ I21 |, [30] ) . 

We discuss two evolutionary strategies. One strategy assumes an overdamped gradient flow 
dynamics and can be viewed as a vanishing viscosity limit (maximum delay convention, e.g. 
[25\ [6]). The other strategy assumes that the system always remains in the global minimum 
of the energy (Maxwell convention). For these two models we establish the existence of the 
thresholds separating the regime of incremental propagation of the decohesion front from the 
regime of a sudden and massive decohesion representing a size effect. When we restrict the 
evolution of the system to the global minimization of the energy, the loading and unloading 
thresholds coincide. When instead we allow the system to follow the maximum delay convention 
and explore the set of marginally stable configurations, the two thresholds become different. The 
comparison with experiments shows that it is the vanishing viscosity solution which reproduces 
the observed adhesion hysteresis most faithfully (p3I Chapter 3]). 

We then develop a continuum analog of our finite dimensional microscopic model, inter- 
preting it as a formal r-limit pQ While the limiting functional, constructed in this way, 
usually captures only the global minima in the original discrete problem, in our case it also 
agrees with a point wise limit and therefore preserves the local minima. To prove this fact we 
perform a systematic study of all metastable solutions of the continuum problem. While the 
continuum model is much more transparent mathematically and allows one to obtain the values 
of all relevant thresholds in explicit form, the strongly discrete limit remains important for some 
applications, in particular, for the modeling of the DNA |22j. 

The primary goal of our simplified model is to elucidate how macroscopic responses de- 
pend on the microscopic parameters. The continuum version of the model contains only one 
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non-dimensional micro-parameter v measuring the propensity of the adhesive layer to dynamic 
snapping. We show that the degree of localization of the decohesion front increases as v decreases 
which is revealed macroscopically as a transition from quasi-ductile to quasi-brittle behavior. In 
the limit v — > 0, which corresponds to the disappearing of an internal length scale, we obtain an 
explicit expression for the critical debonding force, which is in principle a measurable parameter 
[28j . In general, we expect that the achieved parametric control in the simplified microscopic 
setting can be used in the design of the prototypical molecular devices. We are fully aware, 
however, that there is long way between the toy models of the type considered in this paper and 
the realistic description of the particular biological systems (gecko, DNA, etc.) 

The paper is organized as follows. In Section [1] we introduce our discrete model and present 
an analytical description of all stable and metastable configurations corresponding to a given 
value of applied displacement. In Section[2]we derive the T-limit of the discrete model and classify 
the local minimizers of the limiting problem. In Section [3] we study two different responses of the 
discrete model to monotonous and cyclic loading, one overdamped and another nondissipative. 
In the dissipative case we compute the associated heat to work ratio and construct the hysteresis 
loops. Finally, in Section 0] we show how the main features of the cohesion/decohesion hysteresis 
vary as one goes from discrete to continuum description and present the results of the detailed 
parametric study of the model. In the Appendices we collect mathematical results of technical 
nature. 



Consider a discrete chain containing n+1 points which are connected by linear elastic springs 
with reference length / = L/n. Each point is also connected to a rigid substrate by a breakable 
spring. In this maximally simplified setting one can deal with two prototypical problems: pull 
out test (e.g. [17]) and pull off or peeling test (e.g. [SO])- For determinacy, we shall focus on the 
peeling problem and assume that the points move orthogonally to the substrate (see Figj2]). 

Denote by Uj the vertical displacements of the particles from their reference positions. The 
elastic energy of the connecting linear springs can be written as 



where G is the {shear) modulus and <5j = (ui + \ — Ui)/l. For the energy of the breakable springs 
we assume the simplest form 
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FIGURE 2. Schematic representation of the discrete model. The location of the decohesion 
front is marked by the integer valued parameter £. In the horizontal direction the particles are 
separated by the fixed distance I — L/n. 



where E is the longitudinal elastic modulus, Wi = Ui/u r are the normalized displacements, and 
u r is the breaking threshold. The total energy of the chain can be written as 



L n+1 



* = + £<K^)). (i.3) 

71 i=i i=i 

We load the chain in a hard device, meaning that the (normalized) displacement d of the first 
spring is prescribed 

wi = d > 0. (1.4) 

The energy can be rewritten in a more compact form. To this end we introduce a vector 
distinguishing 'bonded' and 'de-bonded' springs 

f if Wi < 1 

X« = (1-5) 
I 1 if Wi > 1, 

and construct the (n + 1) x (n + 1) diagonal matrix 

£> = diag(x(l),...,x(n + l)). 

In what follows it will be convenient to use the following notations: w G M n+1 - the displacement 
vector, 1 6 iR n+1 - the vector with (l)j = 1, and i\ - the first vector of the canonical basis in 
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M n+ . We also introduce the (n + 1) x (n + 1) tri-diagonal matrix 



1 -1 
-1 2 -1 
-1 2 











-1 

-1 2 



-1 

-1 2 -1 
-1 1 



By using the above notations, we can rewrite the dimensionless total energy (11, 3D in the form 

J D (™):=-^ = ^(Bu>. «; + £), (1-6) 

where 

B = I - D + n 2 u 2 A (1.7) 

and ^ is the total number of debonded springs. The dimensionless energy (11.61) depends on two 
scaling parameters: n and 

u r G 

— fVF _ (L8> 

The main physical nondimensional parameter of the problem, z/, implicitly characterizes the 
toughness of the breakable bonds. In particular, by decreasing v we increase the cohesion energy. 
The geometrical parameter n is a measure of discreteness and n — > oo would mean for us the 
'macroscopic' or 'continuum' limit (see Section [2]). 

To find the equilibrium state of the chain at a given d, it is natural to first minimize 
the elastic energy at a fixed configuration of debonded springs D. We obtain the following 
minimization problem 

Mm {j D (w) | w G M n+ \w ■ i x = d) . (1.9) 
The necessary conditions of equilibrium can be written as 

Bw-ani 1 = 0, (1.10) 

where 



EL 



A 



is the Lagrange multiplier, representing the external force A acting on the first point of the chain. 
The stability of these equilibrium configurations is ensured by the positive definiteness of the 
Hessian matrix B which immediately follows from (jl.7p . We can then conclude that all solutions 
of (jl.lOp are local minima of the energy. 

The linear equations (jl.lOp can be solved formally which allows us to represent all metastable 
equilibrium configurations by the formulas 
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-d. 



2n V-B n • % x 



(1.12) 
(1.13) 



We observe that the variables Wi given by (|1.12p decrease as the index i increases. This 
follows from the fact that the column elements in the inversions of diagonally dominant tri- 
diagonal matrices necessarily decrease (see e.g. |19j). Therefore in each metastable configuration 
it is necessarily the first £ springs that are debonded while the remaining n + 1 — £ are bonded. 
This observation allows one to write the following analytical representation for the displacement 
field (see Appendix A for details) 



Wi 



d-{i-l) 



a 



1, 



cosh[(n + 3/2 — i)rj\ a 
sinh[(n + 1 — £)?]] sinh[r//2] 2nv 2 ' 



£ + l,...,n + l. 



Here 



is the stress, 



a 



K 



Kd 



2nv l 



2£ - 1 + coth \ coth[(n + 1 - ^77] 
is the effective elastic modulus, and r\ is one of the two solutions of the equation 

1 



1 + 



cosh [7/]. 



(1.14) 

(1.15) 
(1.16) 

(1.17) 



2n 2 v 2 

Since the equilibrium properties are represented by even functions of 77, the particular choice of 
the solution in (|1.17p is irrelevant. 

The energy of the metastable configurations can be written as 

J = J(d,0 = ^(nKd 2 + 0- (1.18) 

In order to be admissible, the configurations ()1.14|) must respect the compatibility condition 
requiring that all bonded springs have Wi < 1 and all debonded springs have Wi > 1. This 
condition is equivalent to a restriction on £. To obtain this restriction we compute the value of 
the loading parameter d\(£) corresponding to w^+i = 1 and the value g?2(£) corresponding to 
i«£ = 1. We obtain 

2(£-l) 



di(0 



1 + 



1 + 



coth[2] coth[(n - £ + 1)??] + 1 ' 

2^ 

coth[S] coth[(ra - ^ + l)ry] - 1 ' 



(1.19) 
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We call the interval [cii(£), ^2(£)] the stability domain of a solution with a given geometry of a 
crack £. In general, several crack geometries may be compatible with a given load d. The detailed 
study of the obtained solutions, in particular the specialization of the global minimizers, will be 
postponed till Section [3j 

2. Macroscopic problem 

In most applications the parameter n is a large number. Therefore it is of interest to describe 
the continuum limit of the discrete model formulated in the previous section. As a first step we 
can look at the point- wise limits of the discrete solutions (|1.14p as n — ► oo. To compute these 
limits we introduce 

X(i) = (* - l)L/n, 

the coordinate of the ith spring and define the following normalized variables: 

x := X/L, normalized spatial coordinate, 
£ := £/n, fraction of debonded springs. 

By assuming that in the limit n — > oo the parameter £ is finite we obtain from (|1.15[) 

d(C) = 1 + (2.1) 

It is also easy to see that 

lim dt(0 = lim d 2 (() = d((). (2.2) 

n— >oo n— >oo 

This means that for each value of £ the stability domain shrinks in the continuum limit to a 
point. For the limiting displacement field we obtain 



f d — if x € (0, C), 



w^(x) 



V 2 ' 



coshfifl-x)) , , ( 2 - 3 ) 

ifxG(Cl). 



[ cosh(i(l-C)) 

The value of the continuum energy of a metastable state can now be computed from the formula 

J = J(0 = l(t(l + ^)+<r)- (2-4) 
Here we used the limiting relation between the stress and the length of the debonded region 

a = znanh ■ ( 2 - 5 ) 

To interpret these results correctly, we can independently look at the limit of the variational 
problem (II. 9ft as n — ► oo. To this end we can define the space of piecewise constant functions on 
(0,1) 

A n (0,l) =: {5>?V M) i : a? € R, a? = d| 
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and rewrite the discrete energy functional (II. 3p in the form 



Jn(w) 




™($)-<^K 

l/n L 



if we Ax(0, l) 



+00 otherwise in L 2 (0, 1). 



Next we can define A* n as the subset of the functions w G -ff 1 (0, 1) such that there exists w G 
which satisfies 



u/(x) 



^ l/n A [*-i.0±' 

4 = 1 



iw(0) = d. 



Clearly, w( 



w(j:) and we rewrite the original functional in the form 
n+l / / . \ \ . „i 



-y 



i=l 



+ 







-f-w') dx, 
Li 



(2.6) 



(2.7) 



where now w G .4*. It is easy to see that all (local and global) minimizers of J n on A n can be 
described in terms of the corresponding minimizers of J n on .4* which makes the two problems 
equivalent. 

We can now study a point- wise limit of the functional (j2.7|) . It is straightforward to see 
that this finite dimensional variational problem converges as n — > 00 to the infinite dimensional 
problem for the continuum functional 

•1 



J(w) 



1 



(<p(io) + <f>{.-j-v/)) dx if w e ^4(0, 1) 



(2. 



+00 otherwise in L 2 (0, 1) 

which is defined on the space A = {w G fl^O, 1) : w(0) = d}. In Appendix B we prove that the 
point-wise convergence automatically implies T- convergence. In particular, this guarantees that 
the global minimizers of (|2.8p can be viewed as the continuum limits of the global minimizers of 

(H2D- 

The next question concerns the status of the local minimisers of (|2.8p . We say that w £ A 
is a local minimizer of J if there exists 5 > such that for every v G A with ||u; — v\\jji < 5 we 
have J(u>) < «/(^)- In the important case d > 1 we can prove (see Appendix C) that w G A is 
a local minimizer of J if and only if it coincides with a solution of the following system: 



w 



in (0,C) 
io(0) = d; io(C) = 1 



in (C,l) 
w(C) = 1, = 



(2.1 



at C = C) where ^ is a local minimizer of the function J(C) = J( w ()- I n the case d < 1 there 
is only one minimum given by the solution of the Euler-Lagrange equation u 2 w" = w in (0, 1) 
with the boundary conditions w(0) = d and w'(l) = 0. In the special case d = 1 there are 
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two solutions, namely, the solution of the differential equation v 2 w" = w and the homogeneous 
solution w = 1 which means that the detached set can be either empty or coincide with the 
whole interval (0, 1). 

The solution of the linear equations (12. 9p can be computed explicitly. It is easy to show that 
they are given exactly by the formulae (j2.3f) . This means that the n — > oo limit of the metastable 
branches of the discrete model coincides with the metastable solutions of the continuum model. 
Therefore besides ensuring convergence of the global minima our pointwise limit also preserves 
the local minimizers. 

3. Dynamic strategies 

In the previous sections we found for each value of the loading parameter d a variety of 
the accessible metastable configurations. Suppose now that the value of the loading parameter 
is changing quasi-statically. Then the choice of a particular local minimum occupied by the 
system at each value of the loading parameter is controlled by dynamics. Of a particular interest 
are the two evolutionary strategies. The first one represents the vanishing viscosity limit of the 
corresponding viscoelastic (overdamped) problem (e.g. [25]). In this case the system stays in a 
given metastable state till it becomes unstable (maximum delay convention). The second strategy 
imposes that the system is always in the global minimum of the energy (Maxwell convention). 
This behavior can be viewed as the zero temperature limit of the hamiltonian (underdamped) 
dynamics. 

3.1. Viscosity solution. Suppose that the parameter d is monotonically increasing starting 
from the value d = with no debonded springs (point O in Figj3]). The 'virgin' branch becomes 
unstable when the first spring detaches at w\ = 1. According to (|1.19p the decohesion starts at 
d = 1^2(0) = 1. The system then switches to a new metastable branch and we assume that in 
this new branch the only first spring, verifying w\ = 1, detaches whereas all other springs still 
remain in the elastic regime (A — B m. Fig|3]). To check that only one spring breaks one has to 
study the global energy landscape and determine the steepest descending paths (see e.g. [24]). 
We observe that in the continuum limit, according with (12. 2D . a single metastable branch can 
be associated to each d and this indeterminacy is automatically overcome. If the displacement 
is increased further, the debonding continues as the second spring reaches the breaking limit 
at d = ^(1) (C in Fig|3]). This pattern repeats itself as the subsequent springs debond one at 
a time. As we see the system follows a 'pinning-depinning' type of evolution with alternating 
slow elastic stages and sudden transitions between different metastable branches. This behavior, 
with the system switching between the branch with £ debonded springs to the branch with 
£ + 1 debonded springs is possible till <i2(£ + 1) > The numerical solution shown in Fig|3] 

shows that there exists a value of the external load, d = d m d, such that for d > d m d the only 
equilibrium solution is the totally debonded configuration, i.e. £ = n. Thus, when d = d mc [ all 
the remaining elastic springs break simultaneously and the system jumps to the fully debonded 
configuration. In the case of infinite n, we shall be able to find the value of the limiting threshold 
d mc i analytically (see Section [2]). 
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FIGURE 3. Overall force-displacement relation for a system with n = 6 springs and v = 0.3 
Viscosity solution (maximum delay strategy) is indicated by the bold lines. 

3.2. Global minimization. To find the global minimum we have to minimize the energy of 
the metastable equilibrium states with respect to the parameter £. Fig|4] shows by bold lines 
the Maxwell path for the same system as in the previous section. We observe the existence of 
another threshold duax < dmd separating the regime with a progressive debonding from the 
sudden jump to the fully detached configuration. Overall, the resulting stress-strain path is 
analogous to the one in the case of the maximum delay convention. In quantitative terms, the 
Maxwell loading path is lower and the transition to the fully debonded configuration is attained 
at a lower assigned displacement. 




FIGURE 4. Overall force-displacement relation for a system with n = 6 springs and v — 0.3 
Global minimum solution (Maxwell convention) is indicated by the bold lines. Compare with 
FigH 
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3.3. Dissipation. To characterize the dissipation associated with the maximum delay strategy, 
consider, for instance, the path b-c-d-e shown in FigHJ We denote by £ the number of debonded 
springs at the starting equilibrium branch b-c. The subsequent branch d-e will then have £ + 
1 debonded elements. According to the maximum delay convention, the system follows the 
equilibrium branch £ until it becomes unstable at d = cfo© (path b-c). Then it switches to a 
new branch with a smaller energy (jump c-d). To find the energy dissipated during this jump 
event we need to compare the energies (I1.18P corresponding to the two branches £ and £ + 1 at 
a fixed displacement d = d. We can write 

AJ(d,0 := J f+1 (d) - J f (J) = ±-(nd 2 (K i+1 - K() + 1). (3.1) 

The first term in the right hand side represents the difference of the elastic energies (area inside 
the triangle O-P-Q in Figj5]). The second term represents the cohesive energy accumulated by 
the system in the transition between the two states (it does not depend on £). The released 
elastic energy is partially accumulated by the system in the form of cohesion energy and the 
rest is dissipated. The dissipation is zero when the energy difference in ()3. 1 j) vanishes which 
corresponds to the Maxwell path. Instead, along the maximum delay path, represented by 
the points b-c-d in FigHJ the system switches to the new branch in a dissipative way (jump 
c-d) and the dissipated energy AJ(d2 (£)>£) ^ s equal to the area C-c-D-d. In general, according 
to the maximum delay convention the area underneath the stress-strain path, representing the 
external work, can be decomposed into the decohesion energy represented in Fig® by the equal 
triangles of unit area (along the Maxwell path at each switching event the increment of the 
decohesion energy has the same magnitude as the increment of elastic energy), the accumulated 
elastic energy represented by the dark grey and the dissipated energy represented by the light 
areas between the maximum delay path and the Maxwell path. 

3.4. Hysteresis. In Fig® and FigJ7]we illustrate the behavior of the system under cyclic load- 
ing. If the Maxwell convention is operative, there is no hysteresis and the system follows elasti- 
cally the same path for loading and unloading (say, path O-A-B-A-0 in FigJHL). 

If the system is unloaded, after the transition to the fully debonded state (d > dMax) the 
crack heals again at d = dMax, with a sudden transformation of the cohesive energy into elastic 
energy. During such event a finite domain of broken springs reconnects simultaneously (path D-C- 
B in FigJTL). Such snaps are indeed observed in experiments, both for loading and unloading (see 
\13\ Chapter 3] and references therein). Experiments show, however, that the detachment and 
reattachment thresholds can be different with the corresponding systems exhibiting an adhesion 
hysteresis. This suggests that the Maxwell strategy may be less realistic than the maximum 
delay strategy. 

Under the maximum delay convention, if the unloading starts before the system reached 
the fully debonded state (d < d m d), the system shows a limited hysteresis (loop O-A-B-C-D-A 
in Fig®,) which disappears in the continuum limit. If we unload the chain inside this hysteresis 



Under Maxwell convention the transition to the fully debonded state also takes place without dissipation and 
at d — dMax(^Max) the decohesion energy equals the elastic energy and ^Max = nKg d%j ax . 
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FIGURE 5. Decomposition of the external work into the dissipated energy (light grey areas 
C-c-D-d, E-e-F-f), the decohesion energy (middle grey area O-A-B-C-D-E-F-O) and the elastic 
energy (dark gray area o-f-g-o); the global minimum response (dashed bold lines); the maximum 
delay response( continuous bold lines). Parameters are the same as in Fig[3] 




FIGURE 6. Partial cyclic loading for a system with n = 30, and v — 0.3). a) Maxwell conven- 
tion, b) maximum delay convention. 
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loop (from, say, a branch £), the system first deforms elastically until d = d\(£) when the 
last broken spring reconnects and the system jumps back to the branch £ — 1. With further 
unloading the system follows this new equilibrium branch until again at d = — 1) another 
springs reconnects and so on. 

In Figd, we illustrate the behavior of the dissipative system during the complete unloading 
from the fully debonded state. We remark that in contrast to the case of small cycle unloading 
the maximal hysteresis is preserved in the macroscopic limit. 



4. Continuum behavior 

We now turn to the study of the continuum solutions (|2.3|) . Using (|2.1|) one can see that 
there exists a critical value d m d such that for 1 < d < d m d the function J(C) from (|2.4j) has two 
non-degenerate critical points where J'(C) = (one stable and one unstable) while for d > d m d 
there are no such critical points. 

One can also check that for d > 1 the derivative J'(0) = J'(l) = 1 — d 2 < 0, which 
means that the function J(C) behaves as shown in the inserts in FigJHJ Notice also that there 
exists another threshold d max such that for d < d max the global minimum is attained at the 
first of the two critical points, whereas for d max < d < d m d the global minimum is attained at 
the boundary of the domain, £ = 0, describing the totally debonded configuration. Moreover, 
this state remains the only minimizer for the whole interval d > d m d- In Figj9] we show the 
stress-strain and energy-strain diagrams illustrating this behavior of the continuum solutions. 

The critical value of displacement d m d can be obtained from the equation d'(C) = ( see 
FigEJ)) which gives the fraction £ m d of debonded springs at d = d m d- We can write explicitly 

Cmd + c(Cmd) ~ Cmd = 0. 
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FIGURE 8 . Phase diagram for the continuum model. Bold lines represent the maximum delay 
convention; bold-dashed lines represent the global minimization strategy (Maxwell convention). 
Inserts show the structure of the function J(C) in the corresponding intervals. Here v = 1. 





FIGURE 9. Equilibrium force and energy for the continuum analog of the system considered 
in Fig(6] Bold lines represent the maximum delay convention; bold-dashed lines - the global 
minimization strategy (Maxwell convention) ; bold continuous lines correspond to unstable equi- 
libria. 



After solving this equation, we can use (|2.ip to find d m d = d{Q m d)- The displacement d max can be 
obtained by first determining the fraction of debonded springs C, ma x which satisfy J(( ma ,x) = ^(1) 
or 

Cmax o I" ^(Cmax) (1 Cmax) — 0. 

Then, using ([2TT]), one can find d max = d(C, max ). 

The overall comparison of Fig. [7] and Fig. [9] shows that the discrete system has a much 
richer set of metastable states (local minima) than its continuum analog. As n increases, we 
observe two major tendencies: some of the branches of the local minima of the discrete system 
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shrink to points representing the local minima of the continuum system whereas some other 
branches simply disappear. On the contrary, the structure of the global minimum path remains 
basically unaffected as n — > oo. 




FIGURE 10. a) Displacement fields at £ = 0.5 and b) force-displacement diagrams: (a) v = 0.1, 
(b) v = 0.25, (c) v — 10. Dashed lines in b) correspond to unstable equilibrium configurations, 
the horizontal plateaux represents the thermodynamic limit u = 0. 




Ur/L 



FIGURE 11. a) Normalization factor <r(0) as a function of v. b) Debonding force A(0) as a 
function of u r /L. 

Due to the relative simplicity of the continuum model, one can study the dependence of 
the response on the remaining nondimensional parameter u characterizing the toughness of the 
adhesion layer. In Fig llOl, we show three displacement fields corresponding to a given size of 
the crack (£ = 0.5) and to different values of the nondimensional parameter v. We observe that 
the 'localization' of the crack tip predicted by our model increases as v decreases. In Fig llOl, we 
represent the force-displacement diagrams generated by the continuum model at different values 
of v. Here the force is normalized by the debonding threshold, corresponding to £ = 

<t(0) = i^tanh — . 
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As the parameter v decreases we observe an interesting evolution of the stress-strain response 
from an instantaneous (brittle) debonding of the whole chain to a (plasticity type) plateaux, in 
the case of a localized tip. In particular, one can see that the ductility of the system, represented 
by the overall size of the adhesion hysteresis, grows as the parameter v decreases. 

For the continuum system (n = oo) the limit v — > can also be interpreted as a "thermo- 
dynamic" limit because simultaneously u r /L — ► and l/L — ► 0. In this case the normalization 
load vanishes, 

v tanh ► 

v 

however, the limit of the actual critical debonding force remains finite 



A(0) = (L/u r )Ea(Q) = 




In particular, the value of the critical force can be computed explicitly 

A c = VEG. 

In the context of the DNA denaturation, where our equilibrium metastable solutions represent 
domain walls connecting bonded and debonded states, the value A c corresponds to the zero 
temperature unzipping threshold [28j . 
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6. Appendix A 



To invert the tri-diagonal matrix B e in (jl.7p we can use iterative formulas from [19]. We 
first relabel displacements as follows 



u>i, £ = 1, displacements of debonded springs 

Wj + £, j = 1, n + 1 — £, displacements of elastic springs 



(6.1) 



and define the vectors u = (v,i) and v = (vj). Consider the first £ equations (jl.lOp corresponding 
to the debonded part of the chain. After rearrangement, these equations can be rewritten as 



Silt: 



nv 



-1 2 







-1 2 







-1 

-1 2 





ui 




" - + nvd~ 


























nvv\ 



(6.2) 
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where we introduced the £ x £ matrix B\ and added d to both sides of the first equation. 
The parameter v± is the deformation of the first bonded spring. Observe that B\ is a Toeplitz 
tri-diagonal matrix which can be inverted explicitly (see e.g. |114 119]) 



l (i + j-\j -i\)(2[j + 2-\j 
v 



■3 



> 13 » 4n(£ + l) 

Since in the right hand side of (|6.2p only the first and the last elements are different from zero, 



we are interested only in 



( 77 , T\ )' (-^l 1 . 

nf nz^(£ + lj 



ni/(£+ 1) 



Using the first equation of f|6.2j) and (jl.4p we obtain 

2 d — t>i 



The remaining equations give 



a = nu 



a 



£ 



d-(i- 1) — 



(6.3) 
(6.4) 



Similarly, we can reformulate the remaining + £ equations corresponding to the bonded 
part of the chain in the form 



Bov 



2 + -4r 



-l 



o 







1 2 + -^ 



-1 2 + ^ 











V2 











d-(£-l) 







r<7 



where we introduced the (n + 1 — £) x (n + 1 — £) matrix _E?2- Once again we have a tri-diagonal 

Toeplitz matrix and since the diagonal elements satisfy (B2)u > 2, i = 1, rt + 1 — £, we can 
write (see again |llj ) 

_ x _ cosh[(n - £ + 2 - | j - i\)rj\ - cosh[(n - £ + 2 - i - j)rj\ 
[ 2 )tJ 2sinh[r/]sinh[(n-£ + 2)T/] 

The parameter r\ is given by the equation (|1.17p and the results do not depend on the choice of 
one of the two solutions of this equation (indeed the equilibrium solutions are even functions of 
rj). As in the previous case, we need only the first and the last columns of the inverse matrix 

sinh[(n — £ + 2 — i)rj] 



(Bn 



in 



( B 2 )i(n-£+l) 



sinh[(n — £ + 2)rj\ 

sinhfi rj] 
sinh[(n — £ + 2)77] 
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By using these formulas we obtain 

sinh[(n + 2 - g - i)rj\ 1 | 

sinh[(n + 2 - £)t7] ^ 2 n . . 

• \r- 1 ( 6 - 5 ) 

+ sinhKn + 2-eh]^' .= l,..-,Ti + l-e. 

7. Appendix B 

Here we prove that the energy J(w) given by (|2.8p is actually the T-limit (with respect to 
the convergence in L 2 (0, 1)) of J n given by (|2.7|) as n — > 00. 

To prove the T-convergence we proceed in several steps. The first step is to find a tight 
lower bound. The following reasoning is standard (see [5]). 

Proposition 7.1. Assume that w n £ .4* and J n (w n ) < C for every n £ N. TTien to 
subsequences w n — > w weakly in -fT 1 (0, 1) and luei. Moreover if w n — > it; weakly in -ff 1 (0, 1) 
and f/ien 

liminf J n (w n ) > J(w). (7.1) 

rt- -»oc 

Proof. The first assertion is trivial since J n (w n ) < C implies ||ioJ,||.z,2 < C" which together with 
Wn(0) = d yields weak compactness of w n in -ff^O, 1) and that any limit point of u) n belongs to 
A. If in addition w n — > u; weakly in -fT 1 (0, 1) and w£i, then the convexity of yields 



liminf / (f)(-fw' n ) dx > / 4>{—^-w ! )dx. 
n->oo J L Jq L 



Moreover by recalling that w n — ► ui in each L p (0, 1) and that 

V -if ( u/„ ( - 

8=1 V V 

are the Riemann sums of the function <p(w) we get 

thus proving the inequality (|7.ip . □ 
The next step is to prove the existence of a recovery sequence. 

Proposition 7.2. Assume that w € A. Then there exists a sequence w n £ A^ such that w n — > w 
weakly in H 1 and 

Jn{w n ) -» J{W). (7.2) 

Proof. Let w & A and define w n £ .4.* as in (|2.6p . It is readily seen by convexity that J n (w n ) — > 
J(w). ' □ 
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Previous statements prove that J = T— lim J n . The relationship between the minimization 

n— >oo 

problems concerning the functionals J n and J is clarified in the next theorem. 
Theorem 7.3. Let w n G A* n such that 

J n {w n ) ~ inf J n -> 0, (7.3) 

then up to subsequences w n — ► u) weakly in i^ 1 (0, 1) and 

J n (w n ) J(w) = inf J. 

Proof. It is readily seen that (|7.3j) yields J n (w n ) < C for suitable C > and by Proposition 17.11 
we get, up to subsequences, w n — > to weakly in -ff 1 (0, 1) and 

liminf J n (u) n ) > J(w). (7.4) 

n^oo 

Let now i»e A Then by Proposition 17.21 there exists a sequence w n G ^4* such that J n {wn) - » 
J(w) and u> n — > to weakly in fl" . Then either J n (w n ) > J n (w n ) or 

Jn(w n ) - ' m fJn -» 

and so 

J(w) > liminf J n {w n ) > liminf J n (w n ) > <^(^) 
that is J(u)) = min J and inf^* J n — ► min J. 

□ 

8. Appendix C 

We recall that w G A is a local minimizer of J if there exists 5 > such that for every v G A 
with 1 1 if — v\\ H i < 5 we have J(w) < «7(u). 
We first show the following result. 

Theorem 8.1. If w is a local minimizer of J then: 

(1) Ifd< \ then {w > 1} = 0; 

(2) Ifd>\ then either {w > 1} = [0, 1] or {w > 1} = (0, C] with ( G (0, 1). 

Proof. We begin with (2). Let d > 1, then > 1} is a non empty relatively open subset of [0, 1] 
and therefore there exists a countable collection of disjoint open intervals of R, say Ij, j G N 
such that 

{ w >i} = {J(i j n[o, i]). 

Assume by contradiction that for every ( G (0, 1), {w > 1} ^ [0, £), then one of the following 
conditions holds true 

i) {«, > 1} = [0, 1] 
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ii) 3 a G (0, 1) such that (a, 1] C {w > 1} and w(a) = 1 

iii) 3 (3, 7 G (0, 1) such that (^,7) C {w > 1} and to(/3) = ^(7) = 1. 

If ii) holds then let 77 G Cq(0, 1), 17 = in [0, a]: since is a local minimizer we get for every 
e > such that £||r/||#i < 5 

< J(w + £77) - J(u;) = e I (v 2 w'rl + wrjl^^^y) dx + o(e) 

J a 



and by letting e — ► we have 



w'r\ 







that is iy" = in (a, 1). Now, since w(a) = 1 and due to the natural boundary condition 
= 0, we get w = 1 in the whole (a, 1), which is a contradiction. Case iii) follows analogously 
and hence 2) is proven. In order to prove 1) suppose by contradiction that {w > 1} 7^ with 
d < 1. Then either ii) or iii) holds true and a contradiction can be obtained also in this case. □ 

We can now study the relation between the local minimizers of the continuum problem 
(I2.8p and the solution of the linear system (12. 9\\ . 



Theorem 8.2. w G A is a local minimizer of J defined by \2. 8\) if and only if there exists 
C G (0, 1) such that 

w" = in (0, C) 

(8.1) 

w(0) = d; w(C) = 1 



and 



v w" = w in (£, 1) 



w(() = l, w'(l) = 
and C is a local minimizer of J(C) = J(wr). 

Proof. By Theorem 18.11 we have that if w G A is a local minimizer of J then there exists £ such 
that w satisfies (|8.ip and (|8.2|) . Moreover, given 5 > 0, there exists a given small enough e > 0, 
such that w £ , the unique solution of 

w" = in (0, C - e) 
w(0) = d; w(C-e) = 1 



and 



w in (C — e, 1) 
w(C-e) = l, w'(l) = 0, 
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satisfies \\w e — w\\jji < 5. This follows from well known results for elliptic equations with variable 
domains (see [3]). Hence J(w £ ) < J{w) and therefore £ is a local minimizer of the function 

j(0 = jk). 

To prove the inverse statement we have to show that if £ is a local minimizer for J then 
u>£ is a local minimizer for J. Let rj > such that for every |£ — Q < rj, J(C) < J(C) : we may 
choose /3 > such that if v € H 1 ^, 1), v(0) = 0, \\v\\ H i < (3, then + v > 1 in [0, ( - |) and 
ui£ + v < 1 in + 2>1]- Hence 

1 l<~i < ,, , ,.o .\ . 1 /" C " + ^ 



J(^ + u) " (^ 2 |Wf + u '| 2 + 1 ) + - / (^ 2 |-u'| 2 + |u A 1| 2 ) dx+ 



1 

+ 2 



l 

(^ 2 |wj + ?/| 2 + \wf + t>| 2 ) da; 
C+i C 



where u(x) A 1 = u(x) if u(x) < 1 and u(x) A 1 = 1 if u(x) > 1. Here u denotes an absolute 
minimizer of 

l 

n ^o/ (^ 2 |u'| 2 + | it A l\ 2 )dx 



2 



among all u € ff 1 ^ - f, C + |) such that u(C ± §) = iw f (C ± §) + u(< ± !)■ 
Therefore by defining 

«; f (s) + t;(s) in [0, 1] \ [C - |, C + I] 
u(x) otherwise 



w*(x 
we get 



J{w£ + v) > J{w*) > J(w c ). 

An argument very close to that used in the beginning of this Appendix shows that there 
exists a unique (* £ (£ — ^ , £ + |) such that tt > 1 in (£ — |, C*) an d ^ < 1 in ((*, ( + ^). Then, 
taking into account that 



"» l(o, c*) 



eargminj^ 2 ^'! 2 + |u;Al| 2 )dx : u> G F 1 (0,C*),^(C*) = l,^(0) = d| 
and 

«>C*| K . 1} e argminjy (z^ 2 |u/| 2 + |u> A l| 2 )cfe : io € ff 1 ^*, 1), u;(C*) = l| , 



since C is a local minimizer for J and |£ — £*| < ?7/2 , we argue 

J(w c ) = J(C) > J(C) = J(wz), 
thus proving the local minimality of w^. 



□ 
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